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Abstract 



In terms of the Interacting Boson Model, shape invariants for the ground 
state, formed by quadrupole moments up to sixth order, are studied in the 
dynamical symmetry limits and over the whole structural range of the IBM-1. 
The results are related to the effective deformation parameters and their fluc- 
tuations in the geometrical model. New signatures that can distinguish vibra- 
tor and 7-soft rotor structures, and one that is related to shape coexistence, 
are identified. 
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Nuclei are often regarded as drops of nuclear matter as in the geometrical model of Bohr 
and Mottelson. Having a view of nuclei as such geometrical objects leads directly to the 
importance of possible deformations of nuclei. The most important deformation of nuclei 
at low energies is the quadrupole deformation to which we restrict our discussion. These 
quadrupole deformations are of special interest as they enable us to make predictions of 
nuclear properties such as energies or E2 transition strengths of the lowest excited states. 

Conversely one can deduce information about nuclear deformations by observing E2 
transition matrix elements. Indeed from a complete set of E2 matrix elements one can 
calculate model independent moments and higher order moments of the quadrupole operator, 
tensorially coupled to a scalar - the shape invariants. Shape invariants were first introduced 
by Kumar [p]] and Cline [f| in the discussion of a large set of E2 matrix elements obtained 
in Coulomb excitation experiments. Calculating shape invariants in the geometrical model 
shows their connection to the deformation parameters (3 and 7 used by Bohr and Mottelson 
or, to be more precise, to effective values j3 e s and 7 e g and the fluctuations of those. Recently 
Jolos et al. have introduced approximation formulae to the lowest shape invariants in 
the framework of the newly developed Q-phonon scheme These approximations now 

make it possible to determine approximate values of the shape invariants from data by using 
only a few absolute B(E2) values. 

This is a substantial result since the advent of radioactive beams opens up entirely new 
nuclear regions for study but, at the same time, the very low intensities of such beams 
means that data will be sparse and that nuclear structure information must be obtained 
from fewer and simpler-to-obtain data. Hence the importance of the approximations to the 
Q-invariants which allow estimates not only of basic deformation parameters such as f3 and 
7, but of higher moments related to the stiffness of the potential in (3 and 7 and to the 
amount of zero point motion. Such information has seldom if ever been available from any 
nuclear data. With these approximation formulae, they are now accessible from simple data. 

It is therefore important to develop a global view of how these shape invariants behave 
as a function of structure so that they can be effectively used as signatures of structure. It 
is the purpose of this Rapid Communication to map out for the first time the behaviour 
of the five essential invariants, as well as several related quantities, over the full range of 
nuclear structure. To do so we will use the algebraic Interacting Boson Model (IBM) PJT0[| 



to study the behaviour of shape invariants in and between the dynamical symmetry limits of 
the IBM. Formulae will be given to transform the shape invariants into effective deformation 
parameters (3 and 7. The values derived from the algebraic model will be compared to values 
in the appropriate limiting cases of the geometrical model. 

Shape invariants are formed by the isoscalar electric quadrupole operator, which is also 
the E2 transition operator in the Consistent Q Formalism (CQF) ||11|| , 



T{E2) =Q = e B Q IBM , (1) 

where Q IBAI is the quadrupole operator in the IBM 

qIBM = Q x = s +d + d + s + X [d + d] {2) (2) 

and eg is the effective boson charge which is fixed for a given nucleus. We define moments 
up to sixth order of the quadrupole operator in the ground state as 
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Q2= _(ot\(Q-QM) (3) 

?3 = m\lQQQ} {0) \ot)\ (4) 

g4 = (Of\(Q-Q) (Q-QM) (5) 

\(ot\(Q-Q) [QQQ] {0) \o+)\ (6) 

f (otliQQQr [QQQ} m \ot) , (7) 



?5 
<l6 



where a dot denotes a scalar product and [QQQ]^ abbreviates the tensor coupling [Q x 
[Q x Q]( 2 )](°). We should note that q 2 is equal to the total absolute E2 excitation strength 
from the ground state 

q 2 = J2B(E2; Of - 2+) . (8) 

j 

q 2 will be the only quantity in our discussion where an absolute value, namely the effective 
boson charge eg, appears. 

With the moments (0-0) we define the relative dimensionless shape invariants by nor- 
malizing to an appropriate power of q 2 

K n = -^for ne {3,4,5,6}. (9) 

The quantities K n do not depend on the effective boson charge e#. The shape invariants K n 
differ from earlier definitions of Jolos et al. || by normalization constants or tensor coupling. 
In the present definitions no value may become infinite and all shape invariants are exactly 
equal to unity in the limit of the rigid symmetric rotor or, in terms of the IBM, the 577(3) 
limit for any boson number N. 

For the calculation of the shape invariants it is convenient to write the expressions 
for the quadrupole moments (§-0) as sums over E2 matrix elements. Therefore the ten- 
sor properties of the quadrupole operator are taken into account and the unity operator 
1 = J2j,i,M \JiM)(JiM\ is inserted between every pair of quadrupole operators. Using the 
Wigner-Eckert theorem and the unitarity relation of Clebsch Gordan coefficients it is possible 
to write the moments q n as 

g 2 = IK II Q II 2+><2+ || Q || o+) (10) 

i 

?3 = l£<0? II Q II 2+><2+ || Q || 2+><2+ || Q || 0+}| (11) 

q±= E(oniQI|2+)(2niQI|o;)(o+||Q||2+) 

_ <2t || Q || 0+) (12) 
^5 = ^/5 I E (0t II Q II 2+) (n || Q || 2+><2+ || Q || 0+) 

■(0+||Q||2+)(2+||Q||0+)| (13) 
qe= i E (0+|IQI|2 J + )(2+||QI|2+)(2+||g||0+) 

i ,j,k,l,m 

•<0+ || Q || 2+}(2+ || Q || 2+}(2+ || Q || 0+) , (14) 
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involving reduced matrix elements between + and 2 + states only. In general only the 
lowest states contribute to the sums because convergence of the Q-phonon expansion of 
nuclear states is fast ||. Matrix elements between nuclear states that differ by several 
Q-phonons are usually small 0. 

In the model of a quadrupole deformed rotor analytical expressions for E2 matrix ele- 
ments and thus for shape invariants can be obtained. In the rigid rotor the shape invariants 
are functions of the fixed deformation parameters (3 and 7. If we assume a non-rigid defor- 
mation, we can give expressions for the shape invariants as 

(?,ZeR 2X 2 



* = UrJ ^UrJ (15) 

({3 3 cos 3 7 ) _ _ .... 

K * = ^3/2 = COS 3 ^ ( 16 ) 

(P 5 cos 3 7 ) 

^5 = ^ 2 )5/2 ( 18 ) 
(l 6 COS 2 3 7 ) 

^6 = , (19) 

explicitly using expectation values of (3 and cos 37. We can define effective values of the 
deformation parameters /3 e g and 7 e g by Eqs. (|1|,|I6|). The parameter 7 eff = |arccosi^3 is 
given in Table | for the appropriate dynamical symmetry limits of the IBM, where 577(3) 
corresponds to a symmetric rigid rotor, 0(6) to a 7-soft nucleus with maximal triaxiality 
and i7(5) to a vibrator. 

The shape invariants are measures of effective deformation parameters and their fluctu- 
ations. This is made more explicit by defining the following quantities as measures of the 
fluctuations of (3 and cos 37: 

a = TT^Ki =K 4 -1 (20) 
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((3^ 

(/? 6 cos 2 37)-(/3 3 cos3 7 ) ^ 2 ( . 

7 = = K 6 -K 3 (21) 

Using expressions (1^-14,20, ^l|) one can analytically calculate the shape invariants in 
the dynamical symmetry limits of the IBM-1. In this paper we will employ the Extended 
Consistent Q Formalism (ECQF) |TJ of the IBM-1, using the IBM-1 Hamiltonian 



H 



ECQF 



(22) 



with Q x taken from Eq. (0). This simple Hamiltonian contains three parameters (a,(,x)- 
While one parameter (a) sets the absolute energy scale, the wave functions depend only 
on two structural constants ((,x)- F° r a given nucleus the boson number N is fixed. The 
ECQF-Hamiltonian covers the three dynamical symmetry limits as indicated in Fig. |l| We 
note that the structural parameter \ appearing in the shape invariants through the E2 
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transition operator is fixed in the SU(3) limit (x=— y7/2) and in the 0(6) limit (x=0) 
while it is unspecified by the Hamiltonian in the {7(5) limit. Therefore shape invariants in 
the i7(5) limit are functions of the structural parameter x- The analytical expressions for 
the shape invariants and the fluctuations in the {7(5) limit as functions of the boson number 
N and the structural parameter x are 

Af ,5 » = -i= |xl (23) 
<< 5 > = 1 ( 1 - JL. ) (24) 




^ = £C^^ +w _ 2+ _pL) (26 , 

< (5> = 1 T ( 3 * 2 + w - 2 ) • (28) 

For completeness we also give K$ as a function of the boson number N in the 0(6) limit 

K o(6) = 1 (iV-2)(iV-l)(iV + 5)(iV + 6) 
6 ~ 3 [iV(iV + 4)] 2 ' 1 ' 

These results enlarge the well known symmetry triangle for wave functions of the IBM-1 
to a structural ECQF-square for shape invariants and thus for the interpretation of nuclear 
shapes. This fact is illustrated in Fig. |]. The use of a similar rectangular representation 
of the parameter space has also been suggested by Bucurescu et al. [|l3|| . Known typical 
examples for particular points of the ECQF-square are, e.g., 172 Yb for £77(3) , 196 Pt for 
0(6), 116 Cd for {7(5) with x=0 and 152 Sm for large x an d moderate values of ( (see [|TJ]-[16 
and discussion below). 



From Eqs. (^ , |25| , |26| , |28|) we note that the necessary extension of the IBM-1 symmetry 
triangle to the ECQF-square is a finite- N-effect, because the shape invariants of {7(5) wave 
functions converge in the limit N — > oo for any value of x- Table | shows the values of the 
shape invariants and their fluctuations in the dynamical symmetry limits of the IBM-1 for 
an infinite boson number N = oo. Only the quantities given in Eqs. (|23|-p9|) depend on the 
boson number N. 

As we would expect the values of 7 e fr in the symmetry limits are 0° and 30° while the 
effective triaxiality fluctuates in the 0(6) and the {7(5) limits. In the 577(3) rigid rotor and 
0(6) 7- soft limits the (3 deformation is rigid while it fluctuates in the {7(5) vibrator limit. 

We have discussed the shape invariants and fluctuations in the dynamical symmetry 
limits of the IBM using analytical expressions. These values provide useful benchmarks for 
the geometrical interpretation of IBM ground state wave functions. However, the dynamical 
symmetry limits of the IBM and the corresponding geometrical models are idealised, ana- 
lytically solvable limits. More accurate descriptions of the low energy structure of collective 
nuclei can usually be obtained by IBM Hamiltonians outside the dynamical symmetry limits. 
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To gain insight in the structure of actual nuclei the quantities of interest have been calcu- 
lated between the symmetry limits, using the ECQF-Hamiltonian fl22[) . The shape invariants 
and their fluctuations have been calculated gridwise over the whole IBM parameter space for 
N = 10 bosons as functions of the structural parameters ( and x- All calculations have been 
performed by diagonalizing the Hamiltonian numerically using the computer code PHINT 
||17|| . Calculations of the shape invariants have been done by a FORTRAN code (QINVAR) 
which evaluates the PHINT output. 

All quantities behave smoothly and one obtains an impression of how the quantities 
vary outside of the dynamical symmetry limits. Fig. |2] represents the numerical results 
of this work presenting the variation of the most important quadrupole invariants over all 
ranges of structure. The behaviour of the invariants, q2, K^-Kq, between the symmetries is 
interesting. Strong variations towards and for deformed nuclei are typical. The invariant 
K4, which is related to fluctuations in (3 via Eq. (^(]), is one of the few observables that can 
distinguish U(5) from 0(6). This can be useful in newly accessible exotic nuclei since K 4 
can be approximately obtained from the simple expression || 

7 B{E2;A+ ^2+) 
Ki l0B(E2;2t^0t)~ 4 ' ( j 

which involves two observables, easily measured, e.g., by Coulomb excitation experiments. 
The approximation ( |30| ) is valid within about 10% for the ECQF-square and boson numbers 
N > 5. This was numerically checked for the whole ECQF-square and for boson numbers 
N = 5,7, 10, 16. For a detailed analysis an experimental value of K^ pr ' can serve as a 
benchmark for starting points of numerical IBM calculations, which can be optimized to 
reproduce the measured transition strengths. The actual value of K± can then be determined 
from the complete set of calculated E2 transition matrix elements. 

For large N, Kq and 07 are quite different in U(5) and 0(6) which is evident from Table 
|. The bottom right panel of Fig. [| shows er 7 , which gives the fluctuations in 7, gridwise 
over the full structural range. Note, however, that Fig. |^ is calculated for a finite boson 
number which lowers the value of K e and <r 7 in the U(5) and 0(6) limits as seen from 
Eqs. (p6| , |28|p9D . In the SU(3) limit 07 vanishes, which characterizes the SU(3) limit as a 
model for a rigid rotor, also in the 7 degree of freedom. In contrast non-vanishing triaxiality 
fluctuations occur in the U(5) and 0(6) limits, indicating that these limits and the whole 
transitional region between them model 7-soft nuclei. 

Finally, we note that all shape invariants, especially 07, change strongly between SU(3) 
and ?7(5)-like values in an unusual region of the IBM-1 parameter space, namely for moderate 
values of ( and x = ~V7/2. Interestingly, this is just the region appropriate to the nucleus 
152 Sm (C=0.57,x=-v / 7/2) MM- The 

case of Sm is currently under active discussion 
and it seems that it shows a certain degree of shape coexistence between spherical and 
deformed shapes with large effective triaxiality [|l(| . 



Above, we discussed the numerical calculation of the exact shape invariants within the 
sd-lBM-1 parameter space, using the ECQF-Hamiltonian (|22|). One aspect of this work is to 
establish the shape invariants as a convenient link between the geometrical model and any 
other nuclear structure model which is able to calculate E2 transition matrix elements. Here 
we have chosen the algebraical IBM. Our ansatz is alternative to the intrinsic state formalism 



by Ginocchio and Kirson [18| which was used much earlier to link the IBM Hamiltonian to 



the geometrical Bohr Hamiltonian. 
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We note that the effective values of the shape parameters [3 e g and 7 e ff do in general not 
exactly coincide with the minima of a corresponding energy surface for the ground state 
in the deformation parameter plane. However, the shape invariants can easily be used to 
compare the predictions from different nuclear models in a geometrically transparent way. 

In principle, the shape invariants can also be measured directly from extensive nuclear 
structure data, providing a direct test of nuclear structure models. Much more intriguing is 
the common case when only a few key observables, like E2 branching ratios from low-lying 
+ states and 2 + states, are known experimentally and when a phenomenological nuclear 
structure model, like the IBM, can be used to extrapolate the data to a complete set of E2 
transition matrix elements. 

To summarize, we have presented analytic expressions for moments up to sixth order of 
the quadrupole operator in the ground state and we have given definitions for the lowest 
shape invariants up to K§. The shape invariants were calculated analytically in the dy- 
namical symmetry limits of the IBM-1. Formulae were given to derive effective deformation 
parameters and their fluctuations from shape invariants, and thus from IBM-1 calculations. 
A study, using the ECQF-Hamiltonian (0), of the behaviour of the shape invariants over 
a full range of structures has been performed for the first time. It shows the smooth but 
yet widely varying behaviour of the invariants. Thus they can be used to determine the 
properties of nuclei by comparing the calculated invariants to experimentally obtained val- 
ues or to results of fits. Moreover, approximate values of these invariants can be obtained 
experimentally simply from B(E2) values involving just the 2f, 2\ and 4f states, and, for 
K 5 and K®, a B(E2) branching ratio from the appropriate excited + states. 

The invariant K^, as well as the fluctuation <r 7 = Kq — i£f , are of special interest as they 
allow to distinguish between 0(6) and U(5) symmetries which can be difficult otherwise 
B|. Finally, the values of a 7 change most rapidly for IBM-1 Hamiltonians that show shape 



coexistence. 
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TABLES 



TABLE I. Shape invariants and fluctuations, calculated for N = oo in the sd-IBM-l for the 
dynamical symmetry limits using analytical expressions. For comparison corresponding values of 
an effective 7 deformation are given. 



N = 00 


U(5) 


SU(3) 


0(6) 


7eff 


30° 


0° 


30° 


K 3 
K 4 
K 5 
K 6 




7 
5 



21 

25 


1 
1 
1 
1 




1 


1 

3 


(7/3 


2 
5 

21 

25 








1 

3 
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FIGURES 
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FIG. 1. The ECQF-square for transition matrix elements in the IBM-1. U(5) converges to one 
single point for N = 00. 




FIG. 2. a)-e): Shape invariants q<i and K3 - Kq c): The invariant K4 also determines the 
fluctuation erg = K4 — 1. f): The fluctuation o" 7 . All quantities shown in a)-f) are calculated 
gridwise over the ECQF parameter space of the IBM for N = 10 bosons. 
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